From GPE to KPZ: finite temperature dynamical structure factor of the ID Bose gas 
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We study the finite temperature dynamical structure factor S(k, cj) of a ID Bose gas using numer- 
ical simulations of the Gross-Pitaevskii equation appropriate to a weakly interacting system. The 
lineshape of the phonon peaks in S(k,io) has a width oc |fc| 3//2 at low wavevectors. This anomalous 
width arises from resonant three-phonon interactions, and reveals a remarkable connection to the 
Kardar-Parisi-Zhang universality class of dynamical critical phenomena. 



The statistical mechanics of low dimensional fluids, 
both quantum and classical, has long been a source of 
theoretical surprises. To give just two examples: 

1. The long-time tail oc t~ d l 2 in the velocity autocor- 
relation function of a <i-dimensional classical fluid 
invalidates hydrodynamics for d < 2 [IH3] . 

2. The Luttinger liquid description [3] provides a uni- 
versal language for ID quantum liquids, with a 
panoply of phases arising upon perturbation [5] . 

Despite these achievements recent developments in the 
theory of ID quantum liquids away from the low energy 
limit make it clear that our understanding of these sys- 
tems is still rather limited (for a review, see Ref . [5] ) . To 
take a simple example, consider the dynamical structure 
factor S(k,uj) that gives the cross section for inelastic 
scattering from the liquid as a function of momentum hk 
and energy Huj transferred. The Luttinger liquid theory 
predicts that S(k, to) consists only of a pair of delta func- 
tion peaks w = ±c|fc|, with c the velocity of sound, corre- 
sponding to an undamped phonon oscillation. At finite 
k, however, one expects this delta function to broaden 
due to interactions between phonons. Attempts to find 
the resulting lineshape using perturbation theory within 
the Luttinger framework are plagued by divergences [7J, 
whose origin we will describe below. As a result, the pos- 
sibility of capturing the relevant physics within the Lut- 
tinger or hydrodynamic formalism is now viewed with a 
degree of pessimism [8]. 

Almost all of the developments reviewed in Ref. [6] 
pertain to zero temperature. In this work we study the 
dynamical structure factor of a ID Bose gas at finite 
temperature. Aside from being of paramount impor- 
tance in real systems, we will show that finite temper- 
ature brings qualitatively new features that cannot be 
interpreted simply as a smearing of the zero temperature 
lineshape. Using analytical arguments and simulations 
of the Gross-Pitaevskii equation (GPE) appropriate to 
weak interactions and finite temperature T, we find that 
the lineshape of the phonon peak in S(k, lo) has a width 
Tfc oc |/c| 3 / 2 at low k (see Fig. [I]). As well as dominating 
any zero temperature structure (generally oc k 2 ) at low 
wavevectors, the 3/2 power is anomalous relative to the 



-in" 



3 





7 

& 

3 i 






2 
1 






-4 


ID IX i', OlE U.Ub UAL 

| ( W - w,)/fc' V2 


J 


I 


A A 



0T000 0.002 0.004 0.006 0.008 0.010 



UJ 

FIG. 1. Dynamical structure factor S(k,uj) of a ID Bose gas 
described by the Gross-Pitaevskii equation for wavevectors 
k = 2irp/L, p = 64,32, 16 (right to left). L = 5 x 2 13 and 
temperature T = 0.005, with length being measured in units 
of the healing length, and energy in units of the chemical 
potential. Inset: Scaling collapse of the phonon peaks using 
the ansatz Eq. |l]). 

k 2 scaling that follows from linearized hydrodynamics [9] . 

This unusual scaling points to a very rich phenomenol- 
ogy. According to a remarkable recent conjecture [10] . 
the long wavelength dynamics of a classical ID fluid at 
finite temperature is in the Kardar-Parisi-Zhang (KPZ) 
universality class describing interface growth [TTlll3| . 
Specifically, the phonon (Brillouin) peaks in S(k, u>) have 
the scaling form at low wavenumber 

c(±) a \ 1 r ( u±c\k\\ 
Sphononfo") « frfps I Ffc 1 (1) 

where fps(x) is given in Eq. (5.7) of Ref. [13]. The mean- 
ing of Eq. (JTJ is that in a frame moving at the speed of 
sound the density fluctuations moving in the same direc- 
tion behave exactly as the fluctuations of the interface 
slope in the KPZ problem. 

There are very few experiments confirming KPZ scal- 
ing to date (151417] . Our hope is that the results of this 
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Letter will lead to its observation in new systems. For ex- 
ample, the structure factor of ID Bose gases of 87 Rb was 
recently measured using Bragg spectroscopy [18], while 
in Ref. |19j the hydrodynamics of supcrfluid Helium in a 
single nanohole was investigated. In the latter case the 
sound absorption coefficient is presumably more accessi- 
ble than the structure factor, and being oc T k displays 
the same anomalous scaling. 

Hydrodynamic description. Our starting point is the 
classical Hamiltonian describing a ID gas of bosons of 
mass m and interaction parameter g 



H 



d.r 



2m 



9 



(2) 



where the complex field ^(x) obeys the Poisson bracket 
{^(x), ^(y)} = iS(x — y), and we have set H = 1. The 
dynamics of ^(Xjt) is described by the familiar Gross- 
Pitaevskii equation. The mean field description embod- 
ied by the GPE is appropriate when the number of par- 
tides in a healing length £ = (gpom) 1 is large, where 
Po denotes the mean density. This corresponds to 'Lut- 
tinger parameter' K = 3> 1, with c = gp^jm the 
speed of sound in the uniform state. 

After writing the condensate order parameter as 
terms of the canonically conju- 
gate density p(x) and phase 9(x), the Hamiltonian takes 
the form 



H 



dx 



p(d x 6f , {d x ^pf 



2m 



2m 
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2 P 



(3) 



Dynamics near a state of uniform density with p(x, t) — 
Po, 9{x, t) — can be described in the first approximation 
by writing p = p + g, retaining only terms quadratic in 
g and 9 from Eq. ([3| 



Ho = dx 



po {d x ey 

2m 



(d x g) 2 g 2 
8mp 2 



(4) 



H-i is solved by introducing the mode expansions for g{x) 
and 9{x) for a system of length L 



g{x) 



6(x) 




~ Kk (b k e lkx + c.c) 



; Kfc (b k e ikx 



(5) 



After substitution in Eq. Q, e Kfe is chosen to diagonalize 

| 1/2 



2m 



the 



H 2 = Hk E k\b k \ 2 with £ k 

Bogoliubov dispersion relation. At low k £ k — > c\k\ 
0(k 3 ). The deviation from the linear dispersion is due 
to the second term of Eqs. ^ and Q, sometimes called 
the 'quantum pressure'. 



Interactions between the modes are described by the 
anharmonic parts of Eq. ([3]). The most important inter- 
action arises from the first term, and has the form 



H a = dx 



h b k h ( 6 ) 

Kx «2 «3 




^ \/ 32Lpnm 

fei+fc 2 +fc 3 =0 v 



-b, b\ b 1 . - b\ b , b . 



where for simplicity we have assumed the low k limit for 
e Kk . The difficulty associated with a perturbative treat- 
ment of this interaction is now apparent. Substituting 



the time dependence b k — > b k e~" ltk associated with H2, 
we see that the second and third terms of Eq. ^ are res- 
onant for purely linear dispersion when all three modes 
move in the same direction (in quantum mechanical lan- 
guage energy and momentum conservation are simultane- 
ously satisfied). One may object that the 0(k 3 ) deviation 
from linearity at low k due to the quantum pressure term 
removes this difficulty, but the |fc| 3 / 2 broadening that we 
find dominates this effect at low k. The other nonlincar- 
ities arising from Eq. Q are likewise irrelevant in this 
limit. 

The need for a non-perturbative approach was recog- 
nized long ago in Ref. [2U] , where a self-consistent mode- 
coupling (SCMC) treatment of the cubic interaction was 
given, ignoring vertex corrections, and yielded T k cx 
^T\k\ 3 , where T is the temperature (see also Ref. [2"T]). 
The same result was independently rederived much later 
|22) . In contrast, a renormalization group (RG) argu- 
ment based on Galilean invariance (first appearing in the 
related context of the noisy Burgers equation [23]) pre- 
dicts a dynamical critical exponent z = 1 + d/2 for d < 2 
|24j . This suggests that Galilean invariance lies behind 
the success of the SCMC approach, an idea that finds 
support in the analysis of vertex corrections for the case 
of Burgers equation [25] . 

We wish to emphasize that the SCMC theory of 
Rcfs. 20 -22 is an uncontrolled approximation, while the 
RG analysis of Ref. [53] was based on the equations of 
viscous ID hydrodynamics with thermal fluctuations ap- 
pearing as noise sources. It is therefore desirable to study 
the purely Hamiltonian dynamics described by Eq. ([2]). 

Equations of motion. To gain some intuition regarding 
the connection of the GPE to the noisy Burgers equation 
and thence to the KPZ universality class, we discuss the 
equations of motion of the Hamiltonian Eq. ([3]). Ignoring 
the quantum pressure term, these are 



d t p + d x (vp) = 
d t v + vd x v + (g/m)d x p = 



(7) 



where v = d x 9/m is the superfluid velocity. These are 
the continuity and Euler equations for the ID Bose gas, 
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FIG. 2. The characteristic curve X+(t) giving the path of a 
right-moving phonon wavepacket is not straight due to the 
influence of the counterpropoagating waves 



and may be put in the Ricmann form |26) . 

d t {v ± 2c p ) + v ± d x {v ± 2c p ) = 0, 



(8) 



where v± — v±c p , and c p = Cy/pjpo is the speed of sound 
in a frame in which the fluid is locally at rest. Eq. Q tells 
us that v ± 2c p are constant along their respective char- 
acteristic curves X±(t) defined by X±(t) — v±(X±(t),t). 
Alternatively we may write Eq. (fsl) as 



d t v± + v±d x v± = -(d t + v±d x )v : 



(9) 



The interpretation of Eq. (|9| is as follows. Right and 
left moving sound waves propagate through the fluid at 
the velocities v±, but their motion is perturbed by the 
variation in the comoving frame of the velocity of the 
counter-propagating wave. As a result, the characteristic 
curves are not straight lines (see Fig. [2]). We can think of 
Eq. {§) as a pair of driven Burgers equations, in which the 
left moving waves act as a noise term on the propagation 
of the right moving waves, and vice versa. The viscous 
vd 2 v± term is absent for this Hamiltonian system but 
will be generated upon coarse graining. 

It is natural to ask how this situation changes for a 
Fermi gas, which has the hydrodynamic description 



Hi 



Fermi 



i = / dx 



2m 



6m 



(10) 



with the second term representing the Fermi pressure. 
The same analysis now yields the uncoupled Burgers 
equations 



d t v± + v±d x v± = 



(11) 



where v± = v±irp/m are the right and left moving Fermi 
velocities. The characteristics X±(t) are now straight 
lines, and the free Fermi gas therefore represents an ex- 
ceptional fluid in which we expect no anomalous broad- 
ening of the type discussed here. 



Numerical simulations. The GPE is solved using 
the splitting method, whereby 9(x, t) is evolved for a 
timestep r alternately by the kinetic T = \ j dx\d x ^\ 2 
and potential V — \ J dx \ ^\ 4 terms of the Hamiltonian 



T T : #(M) -> e- 4fc2r / 2 *(fc,t) 



(12) 



where ^(k,t) denotes the Fourier transform of ^>(x,t), 
and we now switch to measuring distance in units of the 
healing length £, time in units of (inverse) chemical po- 
tential fx = gp , and \1/ in units of yfpo- Algorithms of this 
type are symplectic. This means that the method exactly 
simulates a Hamiltonian H T with H T — Ha, power series 
in t. The lowest power of r in the series determines the 
order of the method. Two benefits of symplectic integra- 
tors for statistical mechanical simulations are: i) exact 
conservation of phase space volume (i.e. Liouville's the- 
orem is satisfied) and ii) no drift in the energy due to 
exact conservation of H T . 

We use the method V T / 2 • T T • V r / 2 - often called 
'Leapfrog' - which is second order with [27] 



Ht ~ H = 24 (2{T ' {y ' T}} + {y ' {T ' V}}) + ° (r4) 



T 

24 



dx [ P 2 (djdle - (d 2 x ef) - p{d xP f] 

+ 0(r 4 ) 



We display the explicit form of the first correction term in 
H T to demonstrate that the additional terms generated 
by discretizing time are higher order in spatial gradients 
than Eq. and are therefore not expected to change 
the low k behavior. 

The harmonic modes are initially populated accord- 
ing to equipartition, so that the b k are taken as complex 
Gaussian random variables with (|& fe | 2 ) = ^ for tem- 
perature T. This assumes sufficiently weak nonlinearity, 
which is confirmed by the absence of transient behavior 
for the parameter ranges explored, indicating that the 
initial state is close to thermal. 

We choose a spatial discretization scale a = 5 to ensure 
all wavevectors are in the regime of linear dispersion k <C 
The time step of r = 2 is then at the limit of 
stability of the algorithm. Most simulations use systems 
of length L — 2 13 a, but we check that the results are not 
significantly altered for L = 2 14 et (see Fig. |3|. Periodic 
boundary conditions are used throughout. 

At each time step we compute the Fourier components 
of the density p(x, t) 



N-l 

Pk(t) = J2 \^{na,t)\ 2 e 
n=0 



2-irikna 



2tt 



TT 

a 

(13) 
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and to transient rather than equilibrium dynamics, may 
be envisaged. In addition, the challenging problem of de- 
scribing - within a single framework - the finite temper- 
ature phenomena described here, and the zero tempera- 
ture results reviewed in Ref. [BJ, remains to be solved. 
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FIG. 3. Dependence of linewidth on wavevector, for system 
sizes N = 2 13 (red dots) and N — 2 14 (blue crosses), with 
T = 0.005. The data are almost perfectly coincident. The 
blue line is a fit for wavevectors 2irp/ L, p = 2, 4, ... 64 giving 
z = 1.510 ± 0.018. Higher wavevectors show a marked devia- 
tion from 3/2 scaling, and k = 2n / L is omitted because the 
linewidth lies below resolution. 



The resulting time scries is then Fourier transformed to 
give the dynamical structure factor 



S(k,cu) = (\p k:U1 \ 2 ), UJ = 0,± 



2tt 



AT, 



binT 



(14) 



where A^in is the bin size used for the computation of 
the power spectra (at least 2 20 ), and the angular brackets 
denotes an average over ~ 128 runs of length N steps T with 
different random initial conditions. Thus for JVbi n = 2 20 , 
A^teps = 2 21 corresponds to an effective average over ~ 
256 runs, assuming no correlation between the two halves 
of a given run. 

We gather data for wavevectors k — -^2, with p = 
1, 2, 4, . . . 4096. Typical results, displaying good data col- 
lapse (assuming a width scaling as |fc| 3 / 2 ) are shown in 
Fig.[l] For an unbiased test, power spectra are folded as- 
suming symmetry between positive and negative frequen- 
cies, and the phonon peaks fitted to a Lorentzian to ex- 
tract the amplitude, peak frequency, and width. Data are 
shown in Fig. [3] A good fit to the scaling form T k oc \k\ z 
is obtained over 1.5 decades and yields z — 1.510±0.018. 
Significant deviations from scaling are obtained at higher 
wavevectors. 

Conclusion. We have provided convincing analytical 
and experimental evidence of the relationship between 
the finite temperature dynamics of the ID Bose gas and 
the KPZ universality class. Galilean invariance is of 
paramount importance: simulations of wave equations 
with cubic nonlinearity but without Galilean invariance 
show z ~ 1 [25]. Numerous extensions of the results of 
this work to multicomponent (or spinor) quantum fluids, 
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